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Abstract. General theory of elliptic hypergeometric series and 
integrals is outlined. Main attention is paid to the examples obey- 
| ing properties of the "classical" special functions. In particular, an 

■ elliptic analogue of the Gauss hypergeometric function and some 

of its properties are described. Present review is based on author's 
habilitation thesis |Spi7| containing a more detailed account of the 
subject. 
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1. General definition of univariate elliptic 
hypergeometric series and integrals 

Broad definition (n = 1, univariate case) |Spi2[ |Spi4] . 

Formal contour integrals J c A{u)du and series ^ ngZ Cn £1X6 called el- 
liptic hypergeometric integrals and series, if there exist three constants 
lux,uj2,^3 £ C such that 

• A(u + cji) = h(u)A(u), 

where h{u) is an elliptic function of u with periods u>2, UJ3, i.e., 
h(u) is meromorphic and 

h(u + uj 2 ) = h(u + LU3) = h{u), Im(u2/u 3 ) 7^ 0; 

• Cn+i = h(nWi)Cn, 

where h{nuj\) is an elliptic function of n with periods — , — . 

There is a functional freedom in the definition of integrals: A(u) — > 
(p(u)A(u), where <p(u) is an arbitrary cjx-periodic function, if(u + ui) = 
(f(u) (such a freedom is not essential for series). 

Narrow definition of integrals. 

Formal contour integrals J c A(u)du are called elliptic hypergeomet- 
ric integrals, if A(u) is a meromorphic solution of three linear first order 
difference equations 

A(u + u>i) = hi(u)A(u), i = 1,2,3, 

where hi{u) are elliptic functions with the periods uj i+2 (we set 

If all hi{u) 7^ const, then lva.{uji/ujj) 7^ 0, i 7^ j. Interesting situa- 
tions occur when one hi(u) = const, in which case we can have either 
lm(uJi/uJi + i) = or Im^Ui/u^) = 0. For pairwise incommensurate u>i, 
the functional freedom in the definition of A(u) is absent due to the 
non-existence of triply periodic functions. 
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Thus, we have in general three elliptic curves, but only two of them 
are independent. One can consider also elliptic hyper geometric func- 
tions in a more general context, when hi(u) are N x N matrices with 
elliptic function entries. 

It is possible to abandon the requirement of double periodicity of 
h(u) in favor of its double quasiperiodicity similar to the Jacobi theta 
or Weierstrass sigma functions. This leads to a more general family of 
theta hypergeometric series and integrals (theta analogs of the Meijer 
function) |Spi2[ [5pi4] , but we skip their consideration. In the next sec- 
tion we describe certain "classical" special functions of hypergeometric 
type and their elliptic generalizations. 

2. AN OVERVIEW OF CLASSICAL HYPERGEOMETRIC FUNCTIONS 

The Euler's beta integral |AARj 

/ x a -\l - xf-Mx = ^ ^ , Re a, Re(3>0, 
Jo F(a + p) 

determines: i) the measure for Jacobi polynomials; ii) an integral rep- 
resentation for the 2-F1 series. Namely, Jacobi polynomials 

P„ W = %F 1 f-"'" + a + /3 - 1 ;, 
nl \ a 

where (a) n = a(a + l) ■ ■ • (a + n — 1) is the Pochhammer symbol, satisfy 
the orthogonality relations 

(P n ,P m )= [ x a -\l-xf- l P n {x)P m {x)dx 
Jo 

_ S nm T(n + a)T(n + (3) 

~ 2n + a + (3-1 T(n + a + (3 - l)n! ' 
The Gauss hypergeometric function has the form 

a,b \ ^ (a) n (6)n 



x ' n=0 

where we skip for brevity relevant constraints upon the parameters. It 
defines a solution of the hypergeometric equation 

... . ( c a + b — c+l\ . ab . N 

y"{x) +[- + : )y'(x) + — -y(x) = 0, 

\x x — 1 / x(x — 1) 

which is analytical near the origin x = 0. 
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Two integrals described above fit into the general pattern J c A(x)d 



x 



with the kernel A(x) = n^Lil 3 - — x j) aj an d some free parameters Xj 
and ctj. This kernel is characterized by the condition that its logarith- 
mic derivative A'(x)/A(x) = R(x) is a rational function of x. A very 
natural generalization of this criterion consists in the requirement that 
the kernel A(x) satisfies the first order linear finite difference equation 
A(x + uoi) = R(x)A(x) with rational R(x) (such a treatment is valid 
already for the 2 -Pi function via the Mellin-Barnes integral representa- 
tion). By definition, we obtain general plain hypergeometric integrals 
for which 

Ag = TT TV / \ ^ ^ ' <p(x + u h ) = <p{x), 

and r(x) is the Euler's gamma function. 

The Pochhammer series have the form Y2n°=o ° n with 

= R[n) = - — nrr / — : r V, 

c n (n+l)lii=i(« + %) 

which leads automatically to the expression 

E„ fu ,...,U s -l \ ^ (U )n- ■ ■ (U s -i) n n 
c n = s F r [ ; y = > -r/-r — — y . 

n=0 V th,...,v r J f^ j n\{v 1 ) n ---(v r ) n 

These series admit confluence limits like s F r (y) ~ s -iF r {u y) for w — > 
00. Their ^-generalization has by definition c n _|_i /c n = R{q n ) for g G C 
and some rational R(x), which leads uniquely to the series 

fto,---,t s -i \ (*0! Q)n " ' " (ts-li Q)n n 

sVr ;q;y = 2.7 — tt ; 7 > 

where (t; g) n = (1 — t)(l — tq) ■ ■ • (1 — tq"^ 1 ) denotes the g-Pochhammer 
symbol. This definition differs from the one given in [AARl IGR] by the 
inversion q — > q^ 1 and appropriate change of notation for parameters. 
For ti = q Ui , Wi = q Vi and q — > 1~, we formally have s ^Pr(y) — > s F r (y) 
for some renormalized value of the argument y. In a similar way one 
can reconstruct the bilateral series s H r and s ip r . 

Elliptic hypergeometric series, directly derived from the definition 
given in the first section, have the form (the unilateral case) 



r+l 



/ t(), ti, . . . , t r _ \ _ (to)n(tl)n " " " (tr)n n 

\ II)-, ?/!.' ' ' / ' (lilt V. • • • (ill V 



where wq = q (the canonical normalization) and 
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(*)„ = 6(t, tq,..., tq n ~ l -p) = UU W; P) 



8(t;p) = (t;p) 00 (pt 1 ;p) c 



{t;p)c 



nr= (i-*p 



n 



r t- 



w 



!l= o^ »v >, IpI < !; 

the balancing condition. 



The elliptic Pochhammer symbol (t) n (denoted also in some other 
places as (t; q,p) n , 0{t) n) or g) n ) degenerates to (i; g) n for p — ► 0, 

(t)n — * (t] q)n- Therefore for generic fixed tj, Wj we have the termwise 
limiting relation r +iE r —>■ r+i^fr with the balancing restriction indicated 
above (which does not coincide with the balancing condition usually 
accepted for g-hypergeometric series [GR] ) . 



For p 



,,2-Kir 



Im(r) > 0, and any a, u G C, q = e a , we have 



the following relation between 9(t;p) and the Jacobi 6\{x) = 8 1 (x\r) 
function 

oo 

9 1 (au\r) = ~iJ2 ^ k p {2k+1)2/8 q ik+1/2)u 

k=—oo 



Properties 9i(x + 1) 



-0i (a;), 9t(x + r) 



-■KtT — 2'Kix 



6Ax) and 



Q\(— x) = —6i(x) simplify to 9(pz;p) 
For r +\E r series we have 



6(z' 1 ;p) = -z- l 6{z;p). 



Cn+l 



V 



n 

j=0 



6(wjq n ;p) 



h(nu\] 



sua elliptic function of n G C with periods LO%ju)\, uj 3 /uji for 



p 



The integer r + 1 is called the order of h(x) and it counts the number of 
zeros or poles of h(x) inside the fundamental parallelogram of periods. 
Due to the balancing condition, we have an interesting (and useful) 
property 



r+l 



E r 



, t r 



q,p;y 



r+l 



E r 



b > b l ' • • • > b r . -1 



q ,p;y 



to, Hi 

wx,...,w r J \ w x lo- 

in the table below we describe known special functions with prop- 
erties generalizing in a natural way the 2-^1 hypergeometric function 
features. It is rather sketchy and does not pretend on completeness. 
The 2-F1 series is the classical special function investigated by such 
giants as Euler, Gauss, Jacobi, Riemann and many other mathemati- 
cians. Its g-analogue has been proposed by Heine as far back as 1850. 
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However, until the relatively recent time (landmarked by the appear- 
ance of quantum algebras from exactly solvable models of statistical 
mechanics) g-special functions did not attract much attention. 

Chebyshev put forward general theory of orthogonal polynomials 
which played a major role in the search of classical special functions. 
Jacobi polynomials satisfy a three term recurrence relation and the 
hypergeometric equation. The general discrete set of 3.F2 polynomi- 
als was constructed by Chebyshev (I am indebted to R. Askey and A. 
Zhedanov for pointing this fact to me). Their continuous analogues 
and the 302 series generalizations, known as Hahn polynomials, were 
proposed much later. These polynomials satisfy a second order differ- 
ence equation (instead of the differential equation) in their argument 
lying on some non-trivial "grids". The next level of generalization is 
given by the Racah and Wilson polynomials described by special 4F3 
series. In 1985, Askey and Wilson have found [AW] the most general 
set of orthogonal polynomials with the self-duality property. They are 



CLASSICAL SPECIAL FUNCTIONS OF HYPERGEOMETRIC TYPE 



/ Euler \ 
Gauss J 
Jacobi I 

y Riemann/ 



/ Chebyshev \ 

3^ 2 ( 1875 I 



Hahn 



Racah 



rp I 1942 
4^3 I Wilson 



, .-, / Heine 
2^1 I I860 



3^2 (? 9 t) 



/ Askey, \ 



4^3 



1978 

self-dual orthogonal 



I Wilson J 
\ 1985 / 

polynomials 



„T?J Oxgall 
7 1 61 1907 



/ r\ I Jackson \ 
8^1921 ) 



_ / Frenkel, 

-> 10-^9 I Turaev 



summation formulas 



(Wilson \ 
1978 1 
Rahman I 
1986 J 



(Rahman \ 
wSon 
1991 J 



12 



V 1997 



/EE 



2000 

[Sgl4l 

2003 



self-dual biorthogonal rational functions 
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expressed in terms of a special 4^3 series and their argument "lives" on 
the most general admissible grid for polynomials. 

The next level of complexification of functions indicated in the table 
refers to the most general known summation formulas for terminating 
series of hypergeometric type. Sequentially, these are the Dougall's 
7_F 6 and Jackson's sf7 sums going back to the first quarter of the last 
century, and the recent result by Frenkel and Turaev [FT] at the level 
of 10-E9 series to be described below. 

Finally, until very recent time the most general set of known spe- 
cial functions satisfying some orthogonality relations and obeying other 
"classical" properties were given by biorthogonal rational functions re- 
lated to the very well poised $F 8 and 10^9 series. The discrete mea- 
sure functions were discovered by Wilson |Wilj and their continuous 
measure generalizations were derived by Rahman [Rah]. An elliptic 
extension of the Wilson's biorthogonal functions with the key self du- 
ality property was constructed by Zhedanov and the author |SZ1] . The 
Rahman's family of rational functions was lifted to the elliptic level by 
the author |Spi4| . These functions "live" on the grids described by the 
second order elliptic functions — the most general type of grids for ra- 
tional functions admitting a lowering divided difference operator [SZ3J. 
Moreover, in the elliptic case there appeared even more complicated 
objects existing only at this level |Spi4j, which go beyond the space 
of rational functions of some argument and which satisfy unusual two 
index biorthogonality relations. 

There exist also non-self-dual three parameter extension of the last 
row functions described by the very well poised $F S , 10^9, and 12-E11 
series |SZlj . but many of their properties remain unknown. 

In the following we restrict ourselves only to the elliptic hyper- 
geometric functions and for further details concerning plain and q- 
hypergeometric objects we refer to the textbooks |AAR] and |GR| . 
handbook [KS] and the original papers [AWl IRahl IWil] . For a descrip- 
tion of general formal unilateral s E r and bilateral s G r theta hypergeo- 
metric series, see |GRj Spi2 , Spi7|. 



Elliptic hypergeometric integrals are described with the help of the 

bases q,p and 



q = e -'""w^ p 



where q, p, f are modular transforms of q, p, and r. 



Theorem 1. (An elliptic analogue of the Meijer function |Spi4| ) 
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For incommensurate uJi and \p\, \q\, \r\ < 1 general solution of the 
equations 

A(u + (Ji) = hi(u)A(u), i = l,2,3, 

where 



ell. fun-s 


periods 


bases 


moduli 


hi{u) 


UJ 2 , UJ 3 


p 


n = oj 3 /u) 2 


h 2 (u) 


Ll>i, u 3 


r 


T 2 = LOs/uJi 


h 3 {u) 


U)l, U 2 


Q 


r 3 = UJ1/UJ2 



Tl = T 2 T 3 , 



is: 
A(u) 

where 
and 



™ T(t je ^^-p,q) r(t;.e~ 2 -/^;r,g) 
-1-1 T(w j e 27tiu /^]p,q) \\ r(w'e- 2 ™/^i ;r) g) 



o7« 



x constant, 



3 » .'' 3 • - 



rr^n 

3=0 



w 



3 ' 



II II 



T(z;p,q) 



n 

j,k=0 



j=0 j=0 j=0 

1 - z~ 1 p' +1 q k+1 



1 



zpiq k 



q\, \p\ < 1, 



is the standard elliptic gamma function |Rui2j . 
The function T(z;p, q) satisfies equations 

T(z\p,q) = T(z;q,p), T(pq/z;p, q) = 1/T(z;p, q), 
T(qz;p,q) = 6(z; p)T(z; p, q), 
T(pz;p,q) = 6(z;q)T(z;p } q). 

If we denote f(u) = Y{e 2mu ^ 2 ] p, g), then this function solves uniquely 
(up to a multiplicative factor independent on u) the following system 
of three linear first order finite difference equations 



P)f(u); 



f{u + u 2 ) = f(u), 
f(u + uj 3 ) = 6(e 2mu ^;q)f(u). 

There are two choices of parameters with additional nice properties: 

1) 7 = and no t' k ,w' k (the "standard" case \p\, \q\ < 1); 

2) 7 = and m' = m, t'j = rtj, w'j = rvjj (the "unit circle" case). 

In the second case, gamma function factors combine into the modified 
elliptic gamma function introduced in |Spi4| : 



G(u; lo u uj 2 ,uj 3 ) 



Tie 



2mu/u)2 . 



;p,q)T(re 



r, q). 
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This function solves uniquely another system of three equations: 



/(u + wi) = 9(e 2niu/ ^;p)f(u) 



f(u + u 2 ) = e(e 2 ™ u /"i ; r)f(u) 
f(u + uj 3 ) = e- mB ^f(u), 

where 

u 2 u u ui uj 2 1 

B 2 2 [U) = 1 1 1 — . 

tuitu 2 uJi co 2 6uj 2 6u;i 2 

These equations allow us to prove the representation [DS4J 

G(u;uj) = e- niP{u) T(e~ 2m ^;?,p), 



U[U 



2 1 os 



I - , 3uJiUJ 2 UJ3 

related to modular transformations for the standard elliptic gamma 
function |FV] . From this representation it is easy to see that G(u;u>) 
is well defined for \p\, \r\ < 1 and \q\ < 1 (i.e., the \q\ = 1 case is 
permitted in sharp difference from the T(z;p,q) function!). 

Permutations r «-> p and oj\ <-> cj 2 are equivalent. Therefore, we have 

G(u; Ui, u 2 , u 3 ) = G(u; u 2 , uj 1 , u 3 ). 

Due to the property -^(^^=1^ — u) = —P(u), we have the reflection 
equation 

3 

G(a;u>)G(b;uj) = 1, a + b = ^^u k . 

k=l 

In the limit uj 3 — > 00, taken in such a way that simultaneously p,r —>■ 
0, the modified elliptic gamma function is reduced to the "unit circle" 
g-gamma function 



1 , (e 2 ™/^;g) 

hm — = w 1( dj 2 



00 



which remains well defined in the limit \q\ — > 1. This function appeared 
in the modern time mathematics in the work of Shintani [ShiJ as a ra- 
tio of Barnes' double gamma functions [Bar] ; in the works of Faddeev 
and coauthors |Fad[ IFKVj on the modular double of quantum groups 
and quantum Liouville theory; in the work of Jimbo and Miwa [JM] on 
solutions of a g-difference equation and related correlation functions in 
statistical mechanics; in eigenf unctions of the g-Toda chain Hamilton- 
ian [KLSj . In several independent studies it was named as the double 
sign function [Kur] . or hyperbolic gamma function |Rui2l IRui3j . or 
non-compact quantum dilogarithm |FKV] . For the operator algebra 
aspects of this nice function, see |Volj . 
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3. Elliptic functions versus balanced, well poised and 
very well poised hypergeometric functions 

Some convenient terminology. 
Theta functions : holomorphic functions f(x) such that 

f(x + uo 2 ) = e ax+b f(x), f{x + u 3 ) = e cx+d f(x), 

for some a, b,c,d G C. They have a finite number of zeros in the 
parallelogram of periods u>2,u>3, Im^A^) 7^ 0. It is not difficult to 
deduce that 

r 

fix)=e p ^l[e l ix + u J ), Uj eC, 
3=0 

for some polynomial of the second order Pz(x). 

Meromorphic theta functions: ratios of theta functions with different 
parameters r,Uj and P2ix). 

Elliptic functions: balanced meromorphic theta functions 

= tt Oijx + Uj) = -A- 8(tjz;p) 

where p = e 2mT , z = e 2mx , tj = e 2mUj , wj = e 2mVj with the balancing 



constraint 



n$=o h = YYj=o w v or Ei= u 3 = Ei= v i ( mod 1 



guar- 



anteeing that f(x+ 1) = f(x) and f(x + r) = f(x). We can multiply 
these functions by arbitrary independent variable y which is omitted 
for brevity. 

Modular invariant elliptic functions: elliptic functions invariant under 
the action of full PSL(2; Z) group generated by the relations 

f(x; r + 1) = f(x; t), f(x/r; -1/r) = f(x; r). 

Due to the symmetry properties 

6 1 (u\T + l)=e* i '%(u\T), 

- I) = _,(_, r )V2 e ™ 2 /^ i(M | r)) 



elliptic functions are modular if 



E J r =o^ = E; = o^ 2 (mod2r). 



A useful 



form of the second transformation is 

Hip COO • p WQ .UJO+LUQ .U ~U(LOn+LOO:) „ , , 

lif ,e ' = i e m -£^i e m — 44 = e^BiAu) 
B(e "2 ; e -2) 

which indicates that the true modular transformation corresponds to 
the change (u 2 ,w 3 ) — > (-a; 3 ,a; 2 ). 
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Totally elliptic functions: elliptic f(x) which are elliptic also in Uj, Vj 
with the same periods. These are elliptic functions with the con- 
straints Vj = — Uj (modi) or Wj = tj 1 known in the theory of q- 
hypergeometric series as the well poisedness conditions. The balancing 
condition for such well poised elliptic functions YYj=o h = Tij=o ^j 1 * s 
reduced to n^=o*j = i- e - we nave a s ^S n ambiguity! Totally ellip- 
tic functions are invariant under the shifts tj — > ptj (j — 0, 1, . . . , r — 1) 
and z — > pz. Moreover, they are automatically modular invariant and 
satisfy the relation f(—x) = 1/ f(x) (this relation reduces y, the arbi- 
trary multiplier of f(x), to y — ±1). 

We scale now z — > t z and replace parameters tjt by tj (in particular, 
we change t% —> to). As a result, we obtain 



f(z,t) = H 



3=0 



d(tjZ-p) 
B(t~ x z-p) 



n 

3=0 



&(tjz;p) 

e^ot^z-p) 



h(z,t). 



The balancing condition takes now the form YYj=i tj 



±t 



(r-l)/2 




Let 



us take r 
relation 



2k 



n 



2k+l , 



1 odd and resolve the sign ambiguity in favor of the 
Only for this case there are non-trivial sum- 



mation and transformation formulas for series of hypergeometric type. 
In this case h(z, t) is invariant under the shift t — > pt (accompanied 
by the compensating transformation t 2 k+i — > p k t2k+i), i-e. it is an el- 
liptic function of logt with the same periods as for the log z variable. 
Equivalent ly, we have 

f{p 1/2 z,p^H , . . . .A-i.pA) = f(z,t , ...,t r ), 

i.e. there appears interesting symmetry playing with the half period 
shifts. We conclude that the total ellipticity requirement (in appropri- 
ate parametrization) fixes the correct form of the balancing condition 
in the most interesting case of odd r. 

Another important structural constraint leading to interesting ellip- 
tic functions is called the very well poisedness condition. It consists in 
imposing on the well poised elliptic functions of the restrictions 



U-3 



related to the doubling of the Q\(x) function argument. 

We call elliptic hypergeometric series and integrals modular, well 
poised, or very well poised, if the ratios of their kernels c n +i/c n and 
A(u + Ui)/A(u) are modular, well poised, or very well poised elliptic 
functions. It is convenient to introduce special notation for the very 
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well poised elliptic hypergeometric series |Spi3|: 

to, *i, • • • , tr-4, ^v^o, -?v^o ; Qy/to/p, -Qy/P^o 
,qt /h, gtoAr-4, v^o, -\/*o> Vp^o, -y/h/p 



r+1 



q,p; -y 



^ 6{t ;p) 11 



n=0 



m=0 



-(qy) r 



r+l 



v;(* ; *i,---, t r - 4 ;q,p;y), 



m Jn 



where flLi^ 



±tg T 5 ^ 2 g( r 7 )/ 2 (for odd r we assume the positive 



sign, due to the property described above). All known applications of 
these series use a special value of the argument y, y = 1. Therefore, 
we shall drop y in the notation of r +iV r series for this special case. 
For p — > 0, these series reduce to the very well poised r ._iv? r ._ 2 series 
denoted by the symbol r __iW / r _2 in the monograph |GRj . Remarkably, 
the elliptic balancing condition coincides in this case with the usual 
balancing condition accepted for these particular basic hypergeometric 
series [GRl |Spi2l . 

Various forms of the ellipticity requirement provide thus an explana- 
tion of the origin of the notions of balancing and very well poisedness 
for series of hypergeometric type |Spi2j . It is the clarification of these 
points that forced the author to change previous notation for elliptic 
hypergeometric series |Spi2[ |Spi3| . In particular, in this system of con- 
ventions accepted in \GR\ IRos3l |Spi4| , etc the symbol r +iE r used in 
the papers [DSll IKMNOY] ISZ1] should read as r+3 E r+2 or r+3 V r+ 2- 



If we take r 




then 



This is the Frenkel-Turaev summation formula [FT] (for its elementary 
proofs, see, e.g., [Ros3t SZ2J), which is reduced in the limit p — > to 
the Jackson sum for terminating very well poised balanced sf7 series. 

4. The univariate elliptic beta integral 

The elliptic beta integral is the simplest very well poised elliptic 
hypergeometric integral. 

Theorem 2. (The standard elliptic beta integral |Spil|) 
Let ti, . . . , t 6 G C, \tj\ < 1, nj=i tj = P3) an d \p\, \q\ < 1- Then 



K 



Y\l =1 T(t k z;p,q)T(t k z 1 ;p,q)dz 
T(z 2 ;p,q)T(z- 2 ;p,q) z 



Y[ r(t,i fe ;p,g). 

l<?'<fc<6 



where T is the positively oriented unit circle \z\ = 1 and 
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_ (g; 9)00 (p; p)c 

rC — " 

47T« 



The first proof of this integration formula used an elliptic general- 
ization of the Askey's method |Askj which required some contiguous 
relations for the left-hand side expression and Bailey's 2^2 summation 
formula. A very simple proof has been found later on in |Spi6|. 

The elliptic beta integral is the most general univariate beta type 
integral found so far. It serves as a measure in the biorthogonality 
relations for a particular system of functions to be described below. 
After taking the limit p — > 0, our integral is reduced to the Rahman's 
g-beta integral [Rah] 

(<?; q)oo f (z 2 ;q) 00 (z- 2 ;q) 00 (Az;q) 00 (Az~ 1 ;q) 00 dz 



Alxi J-* nSLi(*mz;g)oo(*mZ -1 ;g)oo z 

{tjt k ;q) 00 

where A = Ylm=i ^ m > \^ m \ < 1- This integral determines the measure for 
Rahman's family of continuous biorthogonal rational functions |Rah] . 

If we take now the limit £5 — > 0, then we obtain the celebrated 
Askey- Wilson integral 

(<?; q)oo f (> 2 ;<?)oc0~ 2 ;<?)oo dz {hhht^q)^ 



47li •/Tlll=l( i m^9)oo(^" 1 ;?) D c Z Ill<i<*<4(*J*fc;?)oo' 

determining the measure in orthogonality relations for the most general 
set of classical orthogonal polynomials [AWJ. 

Careful analysis of the structure of residues of the integrand's poles 
allows one to deduce the Frenkel-Turaev summation formula out of the 
elliptic beta integral |DS1] . We suppose that \t m \ < 1, m = 1, . . . , 4, 
\pt 5 \ < 1 < \t 5 \, \pq\ < \A\, A = n s =i*s> an d assume also that the 
arguments of all t s , s = 1, . . . , 5, and p, q are linearly independent over 
Z. We denote C a contour separating sequences of integrand's poles at 
z = t s qip k and A~ x qi +X p +1 , from their reciprocals at z = t~ 1 q~ip~ h , 
Aq~i~ l p~ k ~ l , j, fceN. Then we obtain the following residue formula: 



[ A E {z,t)— = K [ A E { Z ,t)— + Co{t) Y] V n (t) 

Jc z Jt z ~: 



n>0 
I*09"I>1 



with 



n m =i r(* ro « ± \p,q) 
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T(az ± ;p, q) = T(az;p,q)T(az 1 ',p,q), and 

_ (a Ut=i T (tmtf]P,q) 
C ° U T(tf;p,q)T(Atf;p,qy 

» m-,n Q(th 2n ;p) A (U 5 )n 

AJ q 9{%p) U(qt-H 5 ) n 



where we have introduced a new parameter t Q via the relation n m =o tm 
= q. In the limit t 5 t 4 — > q~ N , iV G N, values of the integral on the 
left-hand side of this formula and of the factor co(t) in front of the 
residues sum on the right-hand side blow up, but the integral over the 
unit circle T remains finite. Dividing all the terms by co(t) and taking 
the limit, we obtain the summation formula presented in the end of the 
previous section. 

Using the modified elliptic gamma function it is not difficult to de- 
duce out of the standard elliptic beta integral its "unit circle" analogue 
remaining well defined for \q\ = 1. 

Theorem 3. (The modified elliptic beta integral |DS4j ) 

We suppose that Im{ux/oJ^) > and Im(ujz/uj\) > 0, Im^/a^) > 
and gj G C, j = 1, . . . , 6, Im{gj/uo 3 ) < 0, together with the constraint 

Yf j= i9j = ELi^fc- Then 



where 



/ nf+o S = II G (9j + 9 

(?;?)oo(p;p)oo(r;r)oo 



2(g;g)oo 

Here the integration is taken along the cut with the end points —u^/2 
and UJ3/2. We use also the convention that G(a ± b;u) = G(a + 
b; uo)G(a — b;u). 

If we take Im^) — ► 00 in such a way that p,r — > 0, then this 
integral reduces to a Mellin-Barnes type g-beta integral. More precisely, 
for cg> 12 such that Im(a>i/a>2) > and Re^i/c^) > 0, we substitute 
g Q = J2k=i ^k — A, where A = J2j=i 9j an d a Ppby the inversion formula 
for G(u; uo). Then we set U3 = U0J2, t — > +00, and obtain formally 

S{±2u,A±u;v) du _ 2 (q;q)oo 1^=1 S(A - 9j ; u) 



S(gj + g r , 



where the integration is taken along the line L = iuj^L. Here pa- 
rameters are subject to the constraints Re(gj/uj 2 ) > and Re((^4 — 
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uj\)/uj-2) < 1. This integral was rigorously proven first in [Stoj and a 
quite simple proof was given in |Spi6| in a more general setting. 

5. An elliptic analogue of the 2 F\ function 
We consider the double integral 

nLi F(a j z ± , bjw^ T(cz ± w ± ) dz dw 



l C 2 T(z ±2 ,w ±2 } c 2 Az ± } c 2 Bw ± ) z w ' 

where aj,bj,c G C, A = a^a^, B = b^b^, and C is a contour 
separating converging to zero sequences of poles in z and w from the 
diverging ones, and 

r(*i, ...,t k ) = r(ti;p, q) ■ ■ -V(t k ;p, q). 

Applying the elliptic beta integral formula to integrations with respect 
to z or w (permutation of integrations is allowed since the integrand is 
bounded), we obtain a symmetry transformation for a pair of elliptic 
hypergeometric integrals [Spi4| 

A T(A/ aj ) r nj=irK^,&X) dz 

T{c 2 A/ aj ) Jc T{z ±2 ,cAz ± ,c 2 Bw ± ) z 
= -^j T(B/bj) f H.U r ( a ^' ch i w± ) d z 



r- 



\T(c 2 B/bj) J c T(z ±2 } c 2 Az ± } cBw ± ) z 



This is an elliptic analogue of the four term Bailey transformation for 
non-terminating iov?9 series. It cannot be written yet as some relation 
for infinite uVn elliptic hypergeometric series due to the severe prob- 
lems with their convergence at the boundary values of the argument 

M = i. 

We denote t lj2 ,B = cai,2,3> t 4 = pq/cA, t 5fiJ = fe 1>2 ,3, ^8 = pq/c 2 B and 
introduce the elliptic hypergeometric function — an elliptic analogue of 
the Gauss hypergeometric function 



V(t;,p,q) = k I J r — , [[t j =pq. 



c - • j=1 



Due to the reflection equation for T(z;p, q) function, we have 
V(t]P,q) = I f Titjt^p^), 

t7t8=Pq l<i<k<6 

which is the elliptic beta integration formula (evidently, in this relation 
t 7 and tg can be replaced by any other pair of parameters). 
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In the notation V(t) = V(t;p,q), the transformation derived above 
reads 



v(t)= Yl r(tjt k ,t j+4 t k+4 )v(s) 

l<7<fc<4 



where 



a* = e-Hj (j = 1, 2, 3, 4) _ 7*1*2*3*4 _ r~ro 



s i = £t i U = 5 > 6 > 7 , 8 ) " " V M " V *5*6*7*8 ' 

We repeat this transformation with s 3 , s 4 , s 5 ,s 6 playing the role of 
*i, *2, *3, *4 and permute parameters *3, *4 with * 5 , * 6 in the result. This 
yields 




where T = *i*2*3*4 and [/ = * 5 * 6 * 7 * 8 . 

We equate now the right-hand sides of relations (i) and (ii), express 
tj parameters in terms of Sj and obtain 



m 



V(s)= J] T{s jSk )V(^/pq/s), 
i<i<fc<8 



where y/pq/s = (y/pq/si, VM/sg)- 

Transformations (ii) and (iii) were proven by Rains [Rail] in a strai- 
ghtforward manner using evaluations of determinants of theta functions 
on a dense set of parameters. However, as we just have seen |Spi7| , they 
are mere repetitions of the key transformation (i). 

It is convenient to set temporarily tj = e 2mX3 (pq) l l A . We take vectors 
x G M. 8 and denote as x = ^<=i x « e « their standard decomposition in 
the orthonormal basis e^, (e*, e^-) = Then the balancing condition 
implies Y2^=i x i = which defines a hyperplane Y orthogonal to the 
vector ei + . . . + e 8 . Considering reflections a; — > x — 2(v, x) vj (v, v) 
with respect to the hyperplane normal to some vector v £ Y, it is not 
difficult to see that the transformation of coordinates in (i) corresponds 
to the reflection with respect to the vector v = (Y^i=5 e * — J2t=i e i)/2, 
which has the canonical normalization of the length (v, v) = 2. 

The elliptic hypergeo metric function V(t) appeared for the first time 
in our paper |Spi4| together with the transformation (i). However, 
it was not recognized there that (i) and permutations of parameters 
U <-> tj generate the exceptional E 7 Weyl group of symmetries: the 
function V(t) / rii<fe</<8 \/r(7fc*z) is simply invariant under these trans- 
formations. This fact was understood at the level of series in [KMNOY] 
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(where, actually, only the Eq group is valid since one of the parameters 
is fixed to terminate the series) and for general function V(t) in |Railj . 

For elliptic hypergeometric functions it is convenient to keep two 
systems of notation — the "multiplicative" system, described above, and 
the "additive" one [GR, Spi2, Spi3|. Therefore we define the function 



v{g; u u u 2 , u 3 ) = V^ 2 ™ 91 ^ 2 , e 2ni9a/uJ2 ; e 2 ™ 1 ^ 2 , e 2 ™ 3 ^ 2 ; 



where Y^=i 9k — 2 Y^k=i u k- It will be useful for a description of elliptic 
hypergeometric equation solutions. 



6. Contiguous relations and the elliptic hypergeometric 

equation 

The fundamental addition formula for elliptic theta functions can be 
written in the following form 



a f x y \ a ( x y \ y a x w 

9 xw, —,yz,-;p) -9 xz, -, yw, — ; p) = —9 [ xy, -, wz, —;p 

\ w z / V z w / w \ y z 



where w,x,y,z are arbitrary complex variables. If we denote y = 
t&, w = t 7 , and x = q~Hs, then this identity for theta functions is 
equivalent to the following g-difference equation 



o(t 6 tf- P ) 

9(q-H 8 tf; P y 

t e 9(q-H 8 tt;p) 



A(z, h,..., t 5 , qt 6 , t 7 , q-%) - n ,_J?± __M Z > t) 



t 7 9(q-H 8 t 7 z ;p) 

where A(z,t) = Ylk=i^(^ kZ± ) /^( z±2 ) * s ^ e ^-function integrand. In- 
tegrating now this equality over z along the contour C, we derive the 
first contiguous relation 

t 7 9 (t s t 7 / q, t 8 / qt 7 \ p) V(qt 6 , q~H 8 ) - (t 6 t 7 ) 
= t 7 9 (t 6 t 7 ,t 6 /t 7 ;p)V(t), 

which was used in the first proof of the elliptic beta integral |Spil| . 
Here V(qts, q~ l t$) denotes V(t) with the parameters t e and t 8 replaced 
by qte and q~ l t% respectively and (te <-> t 7 ) means permutation of the 
parameters in the preceding expression. 

In the same way as in the case of series [SZ2j . we can substitute 
symmetry transformation (iii) of the previous section into this equation 
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and obtain the second contiguous relation 

5 

hd(t 7 /qt 8 ;p)Y[6{t 6 t k /q;p)V{q- 1 t 6 ,qt s ) - (t 6 <-> t 7 ) 
k=l 

5 

= t 6 e(t 7 /t 6 ]p)Y[e(t 8 t k ] P )v(t). 

k=l 

An appropriate combination of these two equations yields 

b(t) (U(qt 6 , q-%) - U(tj) + (t 6 <-> t 7 ) + [/(*) = 0, 

where 

m - ^ 



and the potential 

, / x _ 6>(t 6 /gt 8 ,t 6 t 8 ,t 8 /t 6 ;p) -pj 9(t 7 t k /g;p) 
U ~ 9(t 6 /t 7 ,t 7 /qt 6 ,t 6 t 7 /q;p) 11 0(i 8 t fc ;p) 

_ O(qto/t 6 ,t t 6 ,to/t 6 ;p) jj0(q/t 7 t k ;p) 



9(t 6 /t 7 ,qt(i/t 7 ,q/t 6 t 7 ;p) £jj 6(t t k ;p) 

(the second expression is obtained after setting t 8 = p 2 to) is a modular 
invariant elliptic function of variables g±, . . . ,g 7 (tj = e 2 ™ 93 ^ 2 ). 

If we substitute t 6 = az, £7 = a/z and replace [/ (t) by some unknown 
function f(z), then we obtain a g-difference equation of the second 
order called the elliptic hypergeometric equation: 

9(az/qts,atsz,ts/az;p) 



0(z 2 , l/qz 2 ;p) 



j]6(at k /qz;p) (f(qz)-f(z)) 



k=l 

5 



6(a/qt 8 z,at 8 /z,t s z/a;p) , _ x , 

+ — fl(i/zWg;p) M p) (/(? *) - /(*)) 

5 

+ ^(a 2 / g ;p)J]^ 8 ;p)/(z) = 0, 
fc=i 

where t 8 = p 2 q 2 /a 2 nSfc=i*fc- We have found already one functional so- 
lution of this equation U (t) in the restricted region of parameters. The 
second independent solution can be obtained after scaling any of the 
parameters a, t±, . . . , t$ or z by p. We can replace also the standard el- 
liptic gamma functions in the definition of U (t) by the modified elliptic 
gamma functions and get new solutions of the elliptic hypergeometric 
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equation. Indeed, we can rewrite the elliptic hypergeometric equation 

= e 2 ' 



in the "additive" notation tj = e 2 ™ 93 ^ 2 . Then the function 



ui 3 /2 



G(±2x;u) 0J2 



where Ylj=i9i = 2X^=1^; defines its solution linearly independent 
from V(t), provided we impose appropriate restrictions upon the pa- 
rameters. Namely, we should line up sequences of the integrand's poles 
to the left or right of the line passing through the points —003/ 2 and 
o; 3 /2. Evidently, E 7 symmetry remains intact which follows from the 
fact that in the derivation of relevant properties of the V(t) function 
we used only the first (boxed) equation for the elliptic gamma function 
r(z; q,p) which coincides with one of the equations for G(u; o>). Simple 
computations yield the relation 

mod/ . 2u,e 2 ^ p ^Un 9l )) 

V {9 -' u) = l(S-v) (f-f) v ^-^l 

UJ 2{P, P)oo{ r , r )oo 

showing that this solution is proportional to the modular transforma- 
tion of the function v(g; Ui,u 2 , uj 3 ). 

Now we shift gr,s 97,8 + J2k=i an d take the limit Im(a;3) — > 00 
in such a way that p, r — > 0. Then our v mo<i -function is reduced to 

S(±2u, —g 7 ± u, —g$ ± w; ui) du 



s(g;ui,u 2 ) 



H j=1 S(gj±u;uj) u 2 



where Y^=i 9j = 0- This is a g-hypergeometric function which is well 
defined for |g| = 1 and which provides a functional solution of the p = 
degeneration of the elliptic hypergeometric equation. 

It should be noticed that V(t) satisfies not one, but much more equa- 
tions of the derived type due to the permutational symmetry in all its 
parameters, including the equation obtained after the permutation of q 
and p. Most probably there is only one function satisfying all of them, 
since the linearly independent solutions break one of its symmetries, 
E 7 or p <-» q. 

At the level of g-hypergeometric functions, in the limit p — > we 
obtain the equation investigated in detail by Gupta and Masson [GMJ. 
They derived its functional solutions in the form of special combi- 
nations of non-terminating 10^9 series, the integral representation for 
which has been found earlier by Rahman |Rah] and to which our rep- 
resentation for V(t) is reduced in the limit p — > 0. 

In a similar way one can construct contiguous relations for elliptic 
12V11 series with y — 1. Denoting S(t) = 12^11(^0; h, ■ ■ ■ , Q,p), where 
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rim=i t' m = tol 2 anc ^ tm = q n , n £ N, for some m, we have the first 
relation [SZll EZ2] 

0(qk, q\, qt 7 / 1 6 , t 6 t 7 / qt ;p) 



£(t)-£(q~H 6 ,qt 7 



%V*6, q%/t6, to/ h, t 7 / qt ; p) 



5 



nQ(t r ;p) 2 
— — t ;^i, ■ • • ,qt 5 ,t 6 ,qt 7 ), 
r=1 0{qt o /t r ;p) 

and the second one 

0(4-/1- 4-/24- TT+ r C W J 0) ^1) • • • 5 9 J 5, '6) 9 J 7j 

9(t 6 /qt ,t 6 /q 2 t ,te/t 7 ;p) 

, 0(*6;p)nr=lfl(Wg*O;P) P /„2, . - - j. \ 

8(t 7 /qto,t 7 /qHo,t 7 /t 6 ;p) 

9(qt ,qH ;p) U " 

These relations can also be obtained after application of the residue 
calculus similar to the one described above. For this it is necessary 
to take one of the parameters of V(t) outside of the contour C and 
represent this elliptic hypergeometric function as a sum of an integral 
over C and of the residues picked up during this procedure. An accurate 
limit for one of the parameters converting the sum of residues into 
the terminating uVu series brings in the needed contiguous relations, 
which take the described form after changing notation. 

An appropriate combination of these two relations yields 

9(tQ,to/tQ,qto/tQ]p) TT „, , \ (CI 4 -U \ C/4\\ 

I [ 0{qt /t 7 t r ;p) {£ (qt 6 , q t 7 ) - £ (t)) 



9(qt 6 /t 7 ,t 6 /t 7 ;p) 



r=l 



r=l 

5 



9(qt 7 /t 6 ,tr/t a ;p) 

+ e(qt /t 6 t 7 ;p)l[e(t r ;p)£(t) = 0, 

r=l 

which is another form of the elliptic hypergeometric equation. 

7. Applications in mathematical physics 

The theory outlined above did not emerge from scratch. It appeared 
from long time developments in mathematical physics related to classi- 
cal and quantum completely integrable systems. Below we list some of 
the known applications of elliptic hypergeometric series and integrals. 
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(1) Elliptic solutions of the Yang-Baxter equation (elliptic 6j-sym- 
bols) sequentially derived by Baxter |Baxj . Andrews, Baxter 
and Forrester |ABFj . Date, Jimbo, Kuniba, Miwa, and Okado 
|DJKMO] appear to combine into terminating 12 Vn series with 
special discrete values of parameters, as it was shown by Frenkel 
and Turaev in their profound paper [FT] . For a recent work 
in this direction including the algebraic aspects of the ellip- 
tic 6j-symbols, see |Kont IRai2l IRos3l IRos4j . Since solvable 
two-dimensional statistical mechanics models are related to the 
conformal field theory [BM\ IZub] , it is natural to expect that 
elliptic hypergeometric functions will emerge there as well. 

(2) In a joint work with Zhedanov [SZlj , the terminating 12 Vn series 
with arbitrary continuous parameters were discovered as solu- 
tions of the linear problem for some classical integrable system. 
More precisely, these series emerged from self-similar solutions 
of the discrete time chain associated with biorthogonal rational 
functions which generalizes ordinary and relativistic discrete- 
time To da chains. 

(3) As shown by Kajiwara, Masuda, Noumi, Ohta, and Yamada 
[KMNOYJ, Sakai's elliptic Painleve equation [Sakj has a solu- 
tion expressed in terms of the terminating 12V11 series. This 
observation follows from the reduction of corresponding non- 
linear second order finite difference equation to the elliptic hy- 
pergeometric equation. Therefore, V(t) also provides its solu- 
tion. Moreover, the function v(g; u>i, — cu 3 , cu 2 ), well defined in 
the \q\ = 1 region, plays a similar role |Spi7| since it defines 
an independent solution of the elliptic hypergeometric equation 
with the E7 symmetry. More complicated solutions of this equa- 
tion expressed in terms of the multiple elliptic hypergeometric 
integrals were presented by Rains at this workshop |Rai3] . 

(4) Elliptic hypergeometric functions provide particular solutions 
of the finite difference (relativistic) analogues of the elliptic 
Calogero-Sutherland type models [SpTT]. This application is 
outline below and in the last section. 



The original investigations of completely integrable many particles 
systems on the line (or circle) by Calogero, Sutherland and Moser were 
continued by Olshanetsky and Perelomov [OPJ who showed that such 
models are naturally associated with the root systems. Relativistic 
(or finite-difference) generalizations of these models have been discov- 
ered by Ruijsenaars [Ruil] who worked out the A n root system case 
in detail. The corresponding eigenvalue problem is also known to be 
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related to the Macdonald polynomials [Mac]. Inozemtsev |Inoj has in- 
vestigated the most general BC n root system extension of the Heun 
equation absorbing previously derived differential operator models. In 
a further step, van Diejen [Die] has unified Inozemtsev and Ruijsenaars 
models by coming up with even more general integrable model, which 
was investigated in detail by Komori and Hikami [ K H]. A special de- 
generation of this model to the trigonometric level corresponds to the 
Koornwinder polynomials |Kooj . 

The Hamiltonian of the van Diejen model has the form 



n 



+ u(z) 



where u(z) is some complicated explicit combination of theta functions, 
Tjf(. . . , Zj, ...) = /(... , qzj, • • •), and 

A t z \ = Ill=i g Mjiri TT QjtZjZkitZjZ^lp) 
3 ~ 8(z],qz];p) |J- 0( Zj z k , ZjZ^lp) ' 

If we impose the constraint t 2n ~ 2 YtL=i — P 2( fi then the operator 7i 
can be rewritten in the form 

n 

v = E { A ^m - 1) + Mz- 1 )^ - 1)) 

3=1 

up to some additive constant independent on variables Zj (for details, 
see pelimiRuTi] ). 

The standard eigenvalue problem, T>f(z) = Xf(z), in the univariate 
case n — 1 looks like 

I[]=xO{tjz;p) 

-Ff-2 — 2-^ru\<i z ) ~ f(z)) 
0(z 2 , qz 2 ;p) 

Yl^OiUz^-p) 

V[z z ,qz l \p) 

Comparing it with the elliptic hypergeometric equation in the form de- 
rived in |Spi4| , which will be described in the next section, we see that 
they coincide for a restricted choice of parameters t§ — U/q and a spe- 
cial eigenvalue for the Hamiltonian T>, A = — (a similar observation 
has been done by Komori). 
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However, connections between the elliptic hypergeometric functions 
and Calogero- Sutherland type models are deeper than it is just indi- 
cated. Let us introduce the inner product 



where contour C separates sequences of the kernel's poles converging 
to z = from those diverging to infinity. Additionally, we impose 
restrictions upon values of tj and functions ip(z),i()(z), such that we 
can scale the contour C by q and g -1 with respect to the point z = 
without crossing any poles. Under these conditions, the operator V 
formally becomes hermitian with respect to the taken inner product: 
((p,Vip) = (T>ip,i/;). However, this property is not unique — the weight 
function in the inner product can be multiplied by any elliptic function 
p(z), p(qz) = p(z), with an accompanying change of the contour of 
integration. 

In a trivial way, f(z) = 1 is an eigenfunction of T> with the eigen- 
value A = (actually, it solves simultaneously two such equations — the 
second equation is obtained by permutation of q and p). Evidently, the 
norm of this eigenfunction equals to the elliptic hypergeometric func- 
tion, ||1|| 2 = V(t). This relation holds for \p\, \q\ < 1. If we change the 
integration variable in the taken inner product z = e 2mu ^ 2 , then, in- 
stead of f(z) = 1, we could have chosen as an A = eigenfunction of V 
(where the operator T is acting now as a shift, Tv{u) = v{u + Ui)) any 
function h(u) with the property h(u + ui) = h(u), but then the normal- 
ization of this function would not be related to V(t) in a simple way. 
For a special choice of this h(u), we can obtain \\h\\ 2 = v mod (g; to), the 
modified elliptic hypergeometric function for which we can take |g| = 1. 
Equivalently, we could have changed the inner product by replacing the 
standard elliptic gamma functions by their modified version and con- 
sidering the pair of equations Vv{u) = and its w\ <-» uj 2 permuted 
partner. Similar picture holds in the multivariable case considered in 
the end of this paper. 

Because of these relations between V(t) and the Calogero-Sutherland 
type models, it is natural to expect that elliptic hypergeometric func- 
tions will play a major role in the solution of the standard eigenvalue 
problem for the operator D. In particular, we conjecture that the E7 
group of symmetries of V(t) can be lifted to E 8 at the level of uncon- 
strained Hamiltonian TC and that there is some direct relation of this 
model with the elliptic Painleve equation (for this it would be desirable 
to understand an analogue of the Painleve-Calogero correspondence 
principle |LOt [Man] at the level of finite difference equations) . 
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8. BlORTHOGONAL FUNCTIONS 

8.1. Difference equation and three term recurrence relation. 

For n — 0, 1, . . . , we define a sequence of functions |Spi4| 

R n (z;q,p) = 12 V n t s z, —, q ,— ;q,p 

\ £4 tot-4 (- 1 (-4 £2^4 Z 64 

where A = Ylm=o^rn- They solve the elliptic hypergeometric equation 
rewritten in the form 

VJ(z) = 0, V, = V^z)(T - 1) + V^z- 1 )^- 1 - 1) + ^ 

where Tf(z) = f(qz) and 



W) ~ V u ' ^ ' ? ' / d{z 2 ,qz 2 ; P ) ' 



provided we quantize one of the parameters \i = q n ("the spectrum"). 
Equivalently, this equation can be rewritten as a generalized eigenvalue 
problem 

V v f{z) = XDtf{z) 
with the spectral variable lying on the elliptic curve 

where £ and r] are gauge parameters. Out of this representation one 
obtains formal biorthogonality (T n ,R m ) = for n ^ m, where (•, •) is 
some inner product and T n (z;q,p) is a solution of a dual generalized 
eigenvalue problem. 

From the elliptic hypergeometric equation one can derive also the 
three-term recurrence relation 



(y(z) - a n+l )p(Aq n 1 /U)\R n+1 (z;q,p) - R n (z;q,p)^ 

+ (70) - Pn-i) p{q~ n )(Rn-i{z\q,p) - R n (z;q,p)^j 
+ 5(y(z)- 1 (t 3 ))R n (z;q,p) = 0, 
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with the initial conditions = 0, Rq = 1 and 



p(x) 



'' t 4 a:' i 4 x' t *i ' *o*2' *i*2 ' -4 ' A»7 



5 = 6 



0( 

q% q 



qt^x 2 q 2 t&x 2 
~A~l A 

Q 



q_ . *3. 

4 ) , , ) , . ) , , ) ^3*7) j P 
A t l4 t\t{ I2L4 T] 



l(z) = — ; r, <U. 



l{q n /U), /5 n = 7(g n "M). 



9(zr],z/r];p) 

Since the whole z-dependence in this relation is concentrated in the 
j(z) function, R n (z;q,p) are rational functions of j(z) with poles at 
7(2) = ai, . . . ,a n . 

From the general theory of biorthogonal rational functions [Zhel] 
it follows that R n (z;q,p) can be orthogonal to a rational function 
T n (z; q,p) with poles at j(z) = fli, . . . ,(3 n . The involution t 4 — > pq/A 
permutes a n and j3 n , therefore the dual functions are obtained after an 
application of this transformation to R n (z; q,p): 

Ag n-1 



T n (z; q,p) = 12 Vn 



q t ti t 2 z 



U 



q,p 



where the j>dependence in parameters drops out due to the total ellip- 
ticity property (in particular, we have R n {pz; q,p) = R n {z\ q,p)). 

8.2. Two-index biorthogonality. Let us denote the operator T)^ in- 
troduced above as Vjf. Then the product R nm (z) = R n (z;q,p) ■ 
Rm{z;p,q) solves two generalized eigenvalue problems 



m*f{z) = 0, v^f(z) 







q n p m . 



Similar property is valid for the dual 
T n (z;q,p) ■ T m (z;p,q) for a different choice of pa- 



with the spectrum 
product T nm (z) = 
rameters in T>f; p . 

Theorem 4. (Two-index biorthogonality |Spi4| ) 

If we denote 

(g; q)oc(p; p)qq U. 4 m =o r (^, t m z- 1 ) 

Am T(z 2 ,z- 2 ,Az,Az- 1 Y 

n < 



A(z,t) 



ni=o r ( At m 1 ) 

where \q\, \p\ < 1, \t m \ < 1, \pq\ < \A\, then 



dz 

T nl (z)R mk (z)A(z,t) — = h n iAf(t)5 mn 5 k i 



c n 
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where C mn ^i is a contour separating points 



t jP a q b (j = 0, 1, 2, 3), t 4P a - k q b - m ,p a+1 - l q b+1 ~ n /A, a,beN, 

from their z — ► z~ x reciprocals and normalization constants 

h n i = h n (q,p) ■ h m (p,q), 

9(A/qt A ;p)(q,qt 3 /t 4 ,t t 1 ,t t 2 ,t 1 t 2 ,At 3 ) n q~ n 



Only for k — I — there exists the p — > limit and functions 
R n {z; q, 0) and T n (;2; q, 0) coincide with the Rahman's family of contin- 
uous io<^9 biorthogonal rational functions [Rah]. Note also that only 
for fc = / = 0orn = m = 0we have rational functions of some argu- 
ment depending on z\ the general functions R nm (z) and T nm (z) should 
be considered as some meromorphic functions of z with essential sin- 
gularities at z = and z = oo. 

For some quantized values of z and one of the parameters tj the 
functions R n (z;q,p) and T n (z;q,p) are reduced to the finite dimen- 
sional set of biorthogonal rational functions constructed by Zhedanov 
and the author in [SZlJ. They generalize to the elliptic level Wilson's 
family of discrete very well poised w (p 9 biorthogonal functions |Wilj . As 
described by Zhedanov at this workshop [Zhe2] . these functions have 
found nice applications within the general Pade interpolation scheme. 

Functional solutions of the elliptic hypergeometric equation open the 
road to construction of the associated biorthogonal functions following 
the procedure described in [IR] and this is one of the interesting open 
problems for the future. A terminating continued fraction generated 
by the three term recurrence relation described above has been calcu- 
lated in [SZ2]- It is expressed in terms of a terminating 12 Vn series 
and, again, the function V(t) is expected to appear in the descrip- 
tion of non-terminating convergent continued fractions generalizing q- 
hypergeometric examples of [GMj . 

8.3. The unit circle case. In order to describe biorthogonal func- 
tions for which the measure is defined by the modified elliptic beta 
integral, we parametrize tj = e 2nl9j ^ 2 and introduce new notation for 
the functions R n (z; q,p): 



h n (q,p) 



6(Aq 2n /Uq; p) (l/t 3 U, t t 3 , Ms, t 2 t 3 , A/qt 3 , A/qU) 




2ni(w 1 -g 1 -g i )/u} 2 2Tri(u) 1 -g 2 -g4,)/u}2 27ri(^+(n-l)^ 1 -g 4 )/aj 2 




— 2irinu>i/uJ2 p 2-Ki(g:j+u)/u)2 2-Ki(gs— u)/w2 . p 2iriu}\/uJ2 2Ttiuj3/w 2 
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where A = Y^j=a9y Similarly, we redenote the functions T n (z; q,p) as 
s n (u;a>i,a> 2 ,a> 3 ). 

The q <-> p symmetric situation (the standard set of biorthogonal 
functions with \p\, \q\ < 1) is defined as the uj\ <-> kj 3 symmetric product 
of these functions: 

-Rnm(e 27rai/w2 ) = r n (w; c^, uj 2 , u 3 ) • ?m(«; cj 3 , u 2 , ^1), 

with a similar relation for T nm (e 2nlu ^ 2 ) . As described above, we have 
the biorthogonality relations (T nh R mk ) = h n i5 mn 5ku where (1,1) = 1 
coincides with the normalized standard elliptic beta integral with a 
special contour of integration C mn ^i and 

h n i = hn(ui,U 2 ,U3)hl(u}3,U}2,Ul) 

with h n (cui, cu 2 , ^3) = h n (q,p). These functions are modular invariant: 
r n (u; ui, u 2 , uj 3 ) = r n (u; LUi, -u 3 , u 2 ), h n (ui, u 2 , u 3 ) = h n (u>i, -u 3 , u 2 ). 
In the unit circle case we define functions 

r nm( U ) = r n(u; U U U 2 , U 3 ) ■ r m {u] UJ 2 , UJ 1 , L0 3 ), 
S nm( U ) = S n{u] ^1, ^2, ^3) * S m (u\ L0 2 , L0 1 , U 3 ), 

which are now symmetric with respect to the permutations uj 2 <-» <jJ\ 
and n <->• m. These functions satisfy the biorthogonality relations 
(s™° d , r™jf) = h™° d 5 mn 5ki, where (1, 1) = 1 coincides with the normal- 
ized modified elliptic beta integral with the integration contour C mn ,ki 
chosen in an appropriate way and 

h nl° d = M^l,^2,^3) ' ^(^2,^1,^3)- 

In sharp difference from the previous case, the limit p — > (taken in 
such a way that simultaneously r — ► 0, i.e. Ijd.(uj 3 /ui), lm(u 3 /L0 2 ) — > 
+00) exists for all values of indices n, I, k, m and we obtain: 

r nm (u;u 1: u 2 ) = 10 ^9(e 2m(93 - 94)/w2 ;...,e 27r ^- u ^ 2 ;g,g 
x 10 W 9 (e^iaz-ai)/^ . ; e 2ni( g3 -u)/^.~-^ _ 

Their partners from the dual space s nm (u; Ui,u 2 ) are defined in a sim- 
ilar way. These functions r nm (u; ui, uj 2 ) and s nm (u; UJ\, uj 2 ) are not ra- 
tional functions of some particular combination of the variable u for 
n, m 7^ 0. They satisfy the two-index biorthogonality relations 

where v n i are obtained from h™° d after setting Im(a;3/a;i), lva.(uj 3 /uj 2 ) — > 
+00 and (1,1) = 1 coincides with the normalized "unit circle" part- 
ner of the Rahman's integral [Stoj with a special contour of integration. 
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Further simplification of these relations to the Askey- Wilson polynomi- 
als level is highly non-trivial due to some problems with the convergence 
of the integral and requires a thorough investigation. In a similar way 
it is possible to define unit circle partners of the Rains' multivariable 
generalization of the author's univariate biorthogonal functions [Rai2] 
as well as their limiting two-index g-biorthogonal functions. 

9. Multiple elliptic beta integrals 
9.1. General definition. Multiple integrals 

A(«i, ...,«„) dui ■ - -dun, 

where D C C n are some n-dimensional cycles, are called elliptic hy- 
pergeometric integrals if A(u\, . . . ,u n ) are meromorphic functions of 
satisfying the following system of equations 

A(«i, ...,u k + u 1 ,...,u n ) = h {k \u u ...,u n ) A(«i, . . .,u n ), 

where h^(u), k — 1, . . . , n, are elliptic functions of all Uj, i.e., 

h^( Uj + lu 2 ) = h^( Uj + lu 3 ) = hP\u), hn(u 2 /u 3 ) ± 0. 

This is a "broad" definition of the integrals introduced in |Spi4| ; one 
can make it "narrow" by tripling the number of equations for A(u) 
using the shifts by all quasiperiods tu,. 

In order to describe general possible forms of the integrand, we need 
an elliptic extension of the Ore-Sato theorem on the general form of 
terms in plain hypergeometric series (see, e.g., [GGRj ). For all "good" 
known elliptic hypergeometric integrals, the kernels A(w) are equal 
to ratios of elliptic gamma functions T(z;q,p) with an integer power 
dependence on the variables Zj = e 2 ™ 3 ^ 2 . However, in general case 
we can multiply the integrands by elliptic functions of all u/s with the 
periods ui2,u} 3 which do not have such a representation. 

Multiple elliptic hypergeometric series are defined in a similar way. 
It is simply necessary to replace integrals by discrete sums over some 
sublattices of u%, ■■■,u n G Z n keeping other properties of A(w). We 
shall not consider them in the present review. 

The most interesting elliptic hypergeometric integrals are related to 
multiple generalizations of the elliptic beta integral, which are split 
formally into three different groups. Type I integrals contain 2n + 3 free 
parameters and bases p and q and their proofs use in one or another way 
analytical continuation procedure over discrete values of parameters. 
Type II integrals contain less than 2n + 3 free parameters and they can 
be proved by purely algebraic means on the basis of type I integrals. 



/ 
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Finally, type III elliptic beta integrals arise through computations of 
n- dimensional determinants with entries composed of one- dimensional 
integrals. It goes without saying that all these integrals have their 
partners expressed in terms of the modified elliptic gamma function. 

9.2. Integrals for the root system C n . In order to define n-dimen- 
sional type I elliptic beta integral for the root system C n (abbreviated 
as the Cj integral), we take bases \p\, \q\ < 1 and parameters t\, . . . , 

t 2n +4 € C SUch that [1^=1 h = PI anci l*l|> • ■ • > l^2n+4| < 1- 

Theorem 5. (Type I C n elliptic beta integral [DS2J) 



Hi 



C 



11 TV~±2<\ 



n 



dz 



T(zfzt) z 



l<i<j<" V * 3 



n 

l<j<j'<2n+4 



3)i 



where T(z) = T(z; q,p) and 



„c = (p;p)So(g;g)So 



Different complete proofs of this formula were given by Rains [Rail] 
and the author |Spi6| . In the limit p —>■ it is reduced to one of 
the Gustafson results [Guslj . Its modified elliptic gamma function 
partner has been established by the author |Spi6] together with its 
g-degeneration valid for \q\ < 1 (which we skip for brevity). 

Type II integral for this root system (abbreviated as the Cn integral) 
depends on seven parameters t and t m , m = 1, . . . , 6, and bases q,p 
constrained by one relation. It can be derived as a consequence of the 
Cj integral. 



Theorem 6. (Type II C n elliptic beta integral [DM]) 

Let nine complex parameters t,t m (m = l,...,6),p and q be con- 
strained by the conditions \p\, \q\, \t\, \t m \ < 1, andt 2n ~ 1Tlt ' 
Then, 



Um=itm=pq- 



Hi, 



C 



11 r(*±*±\ 11 



l<j<fc<n y 3 k > j=l 



\\[\. II W-X.*. 

j=l \ v ' l<m<s<6 



This is an elliptic analogue of the Selberg integral which appears after 
a number of reductions, the first step being the p — > limit leading 
to one of the Gustafson's integrals [Gus2] . In order to take this limit 
it is necessary to express t$ in terms of other parameters and remove 
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the multipliers pq by the inversion formula for T(z;q,p) (see [DS2J). 
During this procedure the integral takes a less symmetric form — in the 
given form it has the explicit Sq symmetry in parameters (see |Railj ). 
For the modified version of this integration formula valid for |g| < 1 
and its g-degeneration, see [DS4j . 

Presently the author knows only one type III elliptic beta integral 
|Spi4|. It is ascribed to the C n root system (we abbreviate it as the 
Cm integral) and it is computed by evaluation of a determinant of the 
univariate elliptic beta integrals which is reduced to the computation 
of the Warnaar's determinant [WarlJ. We skip it for brevity, but it 
is expected that there are much more such integrals due to the uni- 
versality of the method used for their derivation (see, e.g., |TV] ) and 
existence of several nice exact determinant evaluations for elliptic theta 
functions. 



9.3. Integrals for the root system A n . Classification of the A n el- 
liptic beta integrals follows the same line as in the C n case. We start 
from the description of the simplest type I integral introduced by the 
author in |Spi4] , which we symbolize as Ay . 

Theorem 7. (The A^ integral |5pi4"] ) 



n 



n+1 n+2 



i ^-i^.a n n r(s m zj,t m zj 



l<j<k<n+l^( ZiZ 3 > Z * Z ^ j=lm=l 
n+2 n+2 

= j[T{S8- 1 i Tt- 1 ) J] r( Sk t r , 

fe,m=l 



is dz 

z 



m=l 



where 



Z1Z2 ■ ■ ■ Z n+ \ 



and 



Vn (27ri)"(n+l)! 



with the parameters satisfying the constraints |t m |, \s m \ < 1, m 
1, . . . , n + 2, and I ST = pq, I S = U^=i s m, T = U^=i t m- 



For complete proofs of this formula, see [Rail I Spi6| . Here we have a 
split of 2n + 4 parameters (homogeneous in the C n case) with one con- 
straint into two homogeneous groups with n + 2 entries in each group. 
The p — > limiting value of this integral was derived by Gustafson 
[Guslj . The unit circle analogue together with the appropriate q- 
degeneration valid for |g| < 1 were derived in |Spi6j . Another type 
I A n integral is described below. 
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There are several type II integrals on the A n root system, the first 
of which we abbreviate as Aj) . For its description we define the kernel 

-p/, \ n+l n+1 4 

Ki<i<n+1 l , i ' Z i Z l> j=lk=l i=l 

where t n_1 Ultl h Ylt=i s * = PI and IYj=i z j = 1 - 

Theorem 8. (The A$ integral [Splip 

As a consequence of the Cj and A^p integration formulas, we have 
for odd n 

<\ ^) d i- T ^rf fifing) 

"""" z T ( t 2 A ) fc=iti 

x Yl r(tt 3 t k ) Yl r(tV SiSm ). 

l<j<k<n+l l<i<m<4 

For even n, we have 

„ , n+l 4 

K n / Ag(^ = r(A)nnrN 
n z k=l i=l 

x n r <«*>nSn- 

l<j<fc<n+l i^ 1 ^ ^ 



These formulas contain only n + 5 free parameters. In the p — > 
limit they are reduced to the main result of [GuRj. 

We abbreviate the second type II A n integral as A)/. For its de- 
scription we need the kernel 

p/j. —1 — 1\ n+l 3 

a?»u)= n 'r^.V' ' nn r faw7'>- 

l<i<j<n+l 1 ^ j=l fc=l 

where ten variables p, q, t, s, ti, t 2 , t 3 , si, s 2 , s 3 G C satisfy one constraint 
(ts) n ~ 1 IlLitkS k =pq. 

Theorem 9. (The Af] integral [Spiip 

As a consequence of the Aj , Cj, and Cn integration formulas, we 
have an additional type II elliptic beta integral for the A n root system. 
For odd n, we have 
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A f * (9)/ . dZ , n + l n+1 . T — r n — 1 n — 1 x 

4 Af>(z)—=T(t—,s—) I} Tit—tit^s—SiSk) 

/fn z 

(n+l)/2 3 

x n n r (( ts )^ sfc ) 

j=l j,fc=l 

(n-l)/2 / \ 

x J] r((^') n r^'-V^^V'- 1 ^)) 

j=l \ l<j<fc<3 / 



For even n, we /iai>e 



/ Ag?(2)^ = nr(^, a ^ 



n/2 



X 



xJ]r((^-H A ) n r(t>- 1 s>t i t k ,its>- 1 s i s k ) 

i,k=l l<i<fc<3 



In t/iis and previous theorems we assume constraints on the parameters 
guaranteeing that all sequences of integrands 7 poles converging to zero 
(or their reciprocals) lie within (or outside) ofT. 

This theorem formulas contain only seven free parameters. In the 
p — > limit we obtain one of the integrals in [Gus2] . 

Recently, Warnaar and the author have found a complementary type 
I elliptic beta integral for the A n root system. 

Theorem 10. (The Af ] integral [SWj ) 



n 



YiziZ, 1 , z- 1 z 



l<i<j<n+l 



n n.+3 

nn 



IlUrC^r 1 ) 



1 m =l v m y 1 



II W'T 1 **) 



<?<m<n+3 



where \t k \ < 1 (/c = 1, 2, . . . , n), |s m | < 1 (m = 1, 2, . . . , n + 3), \pq\ < 
\tjS\, S = n^tfi s m , and z 1 --- z n+1 = 1. 

Here we have a split of 2n + 3 independent parameters into two 
groups with n and n + 3 homogeneous entries. This integration formula 
appeared to be new even in the p — > limit as well as in its further 
degeneration to the plain hypergeometric level q — > 1. Its unit circle 
analogue valid for |g| < 1 is constructed in [Spi6|. 
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For each of the described integrals we can apply the residue calculus 
similar to the one described above in the univariate case and derive 
summation formulas for particular multiple elliptic hypergeometric se- 
ries on root systems generalizing the Frenkel-Turaev sum. For the Cj 
integral, the corresponding formula was derived by van Diejen and the 
author [DS3] and its recursive proof was given by Rosengren |Ros2j . 
For the C/j integral, the corresponding sum was conjectured first by 
Warnaar |Warlj . it was deduced from the residue calculus by van Diejen 
and the author [DSlJ and proven recursively by Rosengren [Roslj . The 
Aj resides sum was deduced by the author |Spi4| , leading to an elliptic 
generalization of the Milne's sum [MilJ. Residue calculus for the A^ 
integral performed by Warnaar and the author [S Wj leads to an ellip- 
tic generalization of the Bhatnagar-Schlosser u D n " summation formula 
BS] . These elliptic A^P and A^ summation formulas were proven first 



inductively by Rosengren |Ros2j . A summation formula following from 
the Aj) integral was conjectured by the author |Spi4| , but it still re- 
mains unproven. Residue calculus for the Cm integral is expected to 
lead to a Warnaar's sum | Warl] . but this question was not investigated 
either. 

All the described integrals are expected to serve as measures in the 
orthogonality relations for some biorthogonal functions. A program of 
searching multivariable analogues of the 12 Vn biorthogonal functions 
was put forward in [DS1, Spi2|. The first example of a multivariable 



extension of the author's two-index continuous biorthogonal functions 
was found by Rains |Rai2] on the basis of the C/j elliptic beta integral 
(these functions generalize also the Okounkov's interpolating polyno- 
mials |Oko] ) . 

The notion of root systems provides the main guiding principle in the 
construction of multiple elliptic beta integrals. Although this connec- 
tion is not straightforward, it is natural to expect that there exist other 
such integrals attached, in particular, to the exceptional Lie algebras. 
In this respect it is worth analyzing whether all multiple Askey- Wilson 
type integrals classified by Ito [Ito] admit a further lift up to the levels 
of Rahman's g-beta integral and the author's elliptic beta integral. 

10. Univariate integral Bailey chains 

The Bailey chains techniques is well known as a powerful tool for 
derivation of infinite sequences of identities for series of hypergeometric 
type [AARj . The most general known g-hypergeometric Bailey chain 
was proposed by Andrews [And] . It is related to the Bressoud's ma- 
trix inverse [Ere] and has at the bottom the original constructions 



34 



V. P. SPIRIDONOV 



by Rogers and Bailey used for proving the Roger s-Ramanuj an identi- 
ties |BMj . It was generalized to the elliptic hypergeometric series by 
the author |Spi3| (for some further developments in this direction, see 
[War2j ) . We shall not describe these chains here, although they have 
quite interesting consequences. Instead, we present Bailey chains for 
integrals discovered in |Spi5| . 

Definition. Two functions a(z, t) and (3(z, t) form an elliptic integral 
Bailey pair with respect to the parameter t, if 



dz 

P(w,t)=K I T(tw ± z ± )a(z,t) — . 



Theorem 11. (First integral Bailey lemma |Spi5| ) 

For a given integral Bailey pair a(z,t), (3(z,t) with respect to t, the 
functions 

Y{ts 2 uw~- 



ii \ 1 \^avj i . , 

a (w, st) = ^ U D „ ± , a{w, t), 



s T(t 2 s 2 ,t 2 suw ± ) f T(sw ± x ± ,ux ± ) „, .dx 

P( w > st ) = K Tnnk ±v / v)±2 + 2 2 ±^ 0M — 

T{s 2 ,t z , suw ± ) Jj T(x ±A ,t 2 s 2 ux ± ) x 
where wGT, form a new Bailey pair with respect to the parameter st. 

The proof is quite simple, it is necessary to substitute the key rela- 
tion for (3(x,t) into the definition of (3'(w,st), to change the order of 
integrations, and to apply the elliptic beta integral (under some mild 
restrictions upon parameters). Note that these substitution rules in- 
troduce two new parameters u and s into the Bailey pairs at each step 
of their iterative application. 

Theorem 12. (Second integral Bailey lemma |Spi5|) 

For a given integral Bailey pair a(z,t), f3(z,t) with respect to the 
parameter t, the functions 

.. , r(s 2 t 2 ,M^ ± ) f T(t 2 sux ± , sw ± x ± ) . .dx 
a (w,t) = k — — — / — — a(x,st) — , 

L{S z ,t z ,W ±z ,t z S z UW It )Jj L{SUX ± ) X 

PXw,t) = ^ tU 2 W± l (3(w,st) 
I (ts z uw ± ) 

form a new Bailey pair with respect to t. 

It appears that these two lemmas are related to each other by in- 
version of the integral operator figuring in the definition of integral 
Bailey pairs [SWj . Application of these lemmas is algorithmic: one 
should take the initial a(z, t) and f3(z,t) defined by the elliptic beta 
integral and apply to them described transformations in all possible 
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ways, which yields a binary tree of identities for multiple elliptic hy- 
pergeometric integrals of different dimensions. In particular, the very 
first step yields the key transformation (i) for the elliptic hypergeomet- 
ric function V(t). The residue calculus is supposed to recover elliptic 
Bailey chains for the r +iVr series |Spi3| . We can take the limit p — > 
and reduce all elliptic results to the level of standard g-hypergeometric 
integrals which admit further simplification down to identities gener- 
ated by the plain hypergeometric beta integrals. 

As to the unit circle case, we can start from the relation 



P(v,g) 



W3/2 _ du 

G(g ±v±u;uj) a(u, g) — 
-W3/2 " uj 2 



and apply the modified elliptic beta integral for building needed ana- 
logues of the Bailey lemmas. In this case, the p, r — > limit brings in 
identities for g-hypergeometric integrals defined over the non-compact 
contour L with the kernels well defined for \q\ = 1. 



11. Elliptic Fourier-Bailey type integral 
transformations on root systems 

Similar to the situation with elliptic beta integrals, the univariate 
integral transformation of the previous section has been generalized by 
Warnaar and the author to root systems [SWJ. It appears that in the 
multivariable setting the original space of functions and its image can 
belong to different root systems. 

For the (A, A) pair of root systems, we take the space of meromor- 
phic functions /a(z; t) with A n symmetry in its variables z%, . . . , z n+ i, 
n?=i z j = 1' an d define its image space by setting 



dz 

f A (w;t) = l p&w-^-^fA&t)—, 



D 



where the kernel has the form 



P(z,w;t) 



nSir^ 



> z i z 3, 



In a relatively general situation this map can be inverted explicitly. 

Theorem 13. (The (A, A) transform inversion [SWJ) 

For a suitable n-dimensional cycle D, the inverse of the (A, A) trans- 
form is given by the map 



3(> 
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where it is assumed that functions Ja{x; t) are analytical in a suffi- 
ciently wide annulus encircling T. 

The proof consists in a quite tedious residue calculus with the use of 
the Aj integration formula. 

In the (A, C)-case, we map functions fA{z;t) to its image space be- 
longing to the C n root system: 



dz 



D 



where the kernel has the form 



8 A (z, w;t) 



il 



<i<j< 



n+l~^( Z i Z j i Z i Z ji^ 2 Z{Zj,t 2 Z i Zj ) 



Theorem 14. (The (A, C) transform inversion |SWj ) 

For a suitable n- dimensional cycle D, the inverse of the (A, C) in- 
tegral transform looks as follows 



fA(x;t) 



K 



C 



d c [w, x; t)fc{w; t) — , 



with the kernel 



S c (w,x;t) 



nr=i T ( w t 2 ) iL 



<i<j<n r (^wf) ' 



where it is assumed that functions /a (2,; t) are analytical in a suffi- 
ciently wide annulus containing T. 

Corollary 15. If we choose fc(ML] t) such that the product Sc(w, x; t) ■ 
fc(w] t) is equal to the Cj elliptic beta integral kernel, then the original 
relation fc ~ J d 5a ■ fAdz/z defines the Aj integration formula. 

There are more such Fourier-Bailey type integral transforms with 
explicit inversions some of which still are in the conjectural form. All of 
them can be put into the integral Bailey chains setting yielding many 
infinite sequences of transformations for the elliptic hypergeometric 
integrals on root systems. 



12. Applications to the Calogero- Sutherland type 

models 

After discussing multiple elliptic beta integrals, we would like to 
return to applications of elliptic hypergeometric functions to the Calo- 
gero- Sutherland type models |Spi7| . 
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First, we define the inner product 



dz 

z 

Let us take the Hamiltonian of the van Diejen model [Die] with the 
restriction t 2n ~ 2 Ylm=i t m = P 2 Q 2 



n 

v n = (M*)( T i - 1) + Msr 1 )^ 1 - 1)), 

3 ~ 6(z],qz];p) ~~ @{zj Zk , Zj z^ l ]p) 

Under some relatively mild restrictions upon parameters, this oper- 
ator is formally hermitian with respect to the above inner product, 
((p^iiip) 11 = (Ujj(p,tfj) n , for the weight function 

A n (zf) _ n fr UUm4) 

l<j<k<n 1 Z k > j=i 1 \ Z j > 

Evidently, f(z) — 1 is a A = solution of the standard eigenvalue 
problem T>uf(z) = Xf(z). The norm of this eigenfunction 



i\\ z = v{t-c n ) = K 



is a multivariable analogue of the elliptic hypergeo metric function V(t) 
for the type II C n elliptic beta integral. 

We conjecture that with all multiple elliptic beta integrals one can 
associate Calogero-Sutherland type models in the described fashion. 
Let us take the weight function 

We associate with it the Hamiltonian 

n 

V T = £ (A 3 {z){T 3 - 1) + Mz- 1 )^ - r 
J'=l 

T-r2n+6 /)/, \ n 1 2n+6 

9(z],qz];p) ^9{z j ^p) ^ 

?y 

which is formally hermitian with respect to the taken inner product, 
(<p, Vjip) 1 = (Vjip, ip) 1 , for some mild restrictions upon the parameters. 
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Again, f(z) — 1 is a A = eigenfunction of the operator Vj and its 
normalization 

\\lf = V(t;C I ) = K C n [ A'(z,*)- 

Jin Z 

defines type I generalization of the elliptic hypergeometric function for 
the root system C n . The functions V(t; Cjji) were considered first by 
Rains [Rail] in the context of symmetry transformations for multiple 
elliptic hypergeometric integrals. It is not difficult to define their unit 
circle analogues which also play similar role in the context of Calogero- 
Sutherland type models. 

One can construct analogues of the V(t) function for multiple elliptic 
beta integrals on the A n root system and build corresponding Hamil- 
tonians (all of which coincide in the rank 1 case). Although all these 
models are degenerate — their particles' pairwise coupling constant is 
fixed in one or another way, it would be interesting to clarify whether 
these models define new completely integrable quantum systems. 
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